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From a series of total energy calculations by the full-potential linear muffin-tin orbital method, 
the total energy hypersurface as function of atomic displacements from equilibrium positions has 
been fitted for different F phonon modes in orthorhombic KNbOa. Frequencies and eigenvectors 
of all TO F phonons have been calculated in the harmonic approximation, and in the quantum 
oscillator scheme - for A2 and B2 modes. The microscopic structure of the orthorhombic phase 
^^ ' has been analyzed in a series of supercell calculations for different patterns of Nb displacements, 

Q>^ ' providing indications in favour of the chain structure, with oppositely directed neighboring chains. 
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Z ' I. INTRODUCTION 

^N ' Of three ferroelectric phases known for potassium niobate, the orthorhombic phase is probably the best experi- 

, • mentally studied. This is a room-temperature phase which exists in the temperature range —10 to 225 ''C and is 
►^ \ therefore the most important for any practical applications, e.g., related to the photorefractive properties of KNbOs. 
,_^ . Considerable amount of information has been accumulated from Raman scattering [0,0,p| and neutron scattering B 
,-H ' experiments; local environment of Nb atoms has been probed by EXAFS measurements [g[ ; recent experiments on im- 
T—i , pulsive simulated (time-resolved) Raman scattering ||M| provide additional data on oscillatatory and relaxation-type 
^~~^ ' phonon modes. However, the understanding, on the theoretical level, of fine mechanisms behind ferroelectric atomic 
displacements and phase transitions in KNbOs is still not sufficient. Ah initio calculations which may clarify this 
point have not been up to now performed, to our best knowledge, for the orthorhombic phase, because of the lower 
symmetry and correspondingly higher calculational effort demanded, as compared to the cubic (non-polar) phase and 
Trt ■ the tetragonal ferroelectric phase. 

r^ ' In the present paper, we extend our previous research on the equilibrium geometry of the ferroelectric ground state 

i' , structure of KNbOs ||] and on phonon properties in cubic and tetragonal phases [^ onto the orthorhombic phase of 
'^ ■ this material. The band structure and total energy calculations, based on which the equilibrium geometry and phonon 
^ , prop erties are investigated, have been performed using the full-potential linear muffin-tin orbital code by Methfessel 
O 1^, 11|. The technical aspects of the calculations for KNbOa are discussed in some detail in Ref. |^. 

., Apart from lower symmetry, another difficulty in performing the calculations for the orthorhombic phase is the 

^ much smaller magnitude of relevant energy differences, e.g., related to the displacements within the soft mode, as 
compared to cubic or tetragonal phases. This demands an accuracy of total energy calculation of at least 0.01 mRy. 
In the computer code we applied, that is technically possible due to the use of the Harris energy which coincides with 
cd the Kohn-Sham energy in the first order with respect to variations around the self-consistent density and converges 
much more smoothly than the latter. (More details on the properties of the Harris functional and its applications may 
be found in Ref. |12|]). Also, the convergence of the total energy in the number of special points in the Brillouin zone 
has to be particularly strictly controlled. The corresponding error was always within 0.01 mRy in our calculations 
even for different supercell geometries (with different patterns of the k-points for integration) if they actually refer to 



0\ 



^ 



the same crystal (super) structure |13|. 

The paper is organized as follows: In section ||, the geometry optimization in the orthorhombic phase is discussed 



and compared with the experimental data. In section HI , the results of the frozen phonon calculations in the harmonic 
approximation are presented, and in section ^ - the analysis of phonon frequencies in terms of the quantum oscillator 
model. Finally, in section M different models for the microscopic structure of the orthorhombic phase are considered, 
and the relevance of the chain structures is emphasized. 

II. OPTIMIZED STRUCTURE OF THE ORTHORHOMBIC PHASE 

In discussing the geometry of the orthorhombic phase, we follow the notations introduced by Hewat [IJ]: three 
lattice vectors are a=3.973 A along x = [100] of the cubic aristotypc; 6=5.695 A along y = [Oil] and c=5.721 A 



along z = [Oil]. The positions of atoms and the magnitudes of their off-center displacements in terms of these lattice 
vectors are given in Table |. It would be in principle possible in a first-principle calculation to minimize the total 
energy with respect to all seven independent structure parameters, i.e., three lattice vectors and four internal atomic 
displacements. The optimization of such multidimensional function is however computationally very demanding and 
hardly likely to produce any new essential information. Based on our previous experience, gained from the study of 
tetragonal KNbOa, that the lattice strain found from the total energy minimization is in very good agreement with 
experiment S, and in order to keep the number of independent parameters in the structure optimization manageable, 
we fixed the lattice volume and strain in our frozen phonon calculations and used for the lattice parameters the 
experimental data of Ref. H]. However, all off-center displacements of atoms, which are allowed by symmetry, have 
to be optimized in the calculation, in order to guarantee that for each phonon, the force constants are calculated 
near the total energy minimum. The optimization has been performed by the polynomial fit of the total energy as 
function of four variables; the values of parameters corresponding to the total energy minimum are also given in Table 
fl Similarly to the choice of Ref. [0|, we kept Nb in the center of the cube and allowed all other atoms to relax 
with respect to it. Comparing the calculated and experimental values of A's in Table g, one can conclude that the 
z-displacement of oxygen atoms is fairly well reproduced, including the fact that Oi is slightly farther shifted than 
O2, O3. The y-displacement of O2 is larger than the experimental estimate; it should be noted however that the 
magnitude of this displacement is very small, and of primary importance is the fact that the calculation correctly 
reproduces the sign of this parameter. Finally, the z-displacement of K is overestimated by almost a factor of 3, that 
is identical to what we have found earlier when optimizing the geometry of the tetragonal phase |^. In other words, in 
our calculations K tends to remain coupled to the oxygen sublattice, so that the ferroelectric distortion arises roughly 
due to the pure displacement of Nb. 

III. FROZEN PHONONS 

Combined atomic displacements in the orthorhombic phase which has the space group Amm2 (No. 38 of the 
International tables) are split by symmetry into four one-dimensional irreducible representations. After projecting 
out three translational modes similarly to how it was discussed in some detail by Cohen and Krakauer for cubic 
BaTiOa [15|, one stays with 12 symmetry coordinates St. The choice of St (arbitrary up to constructing other 
linear combinations within each irreducible representation) we used in our total energy fitting and frozen phonon 
calculations, in terms of cartesian displacements of atoms 1 to 5 according to Table ffl, is given in Tabic O, 

We calculated the total energy within each irreducible representation for a sufficient number of different displacement 
patterns as to make a fit by the fourth-order polynomial. Apart from extracting second derivatives for the force 



constants, this fit made possible the quantum oscillator analysis described below in section EV and the optimization of 



finite displacements within the soft mode. A list of polynomial coefficients for the B2 representation, that is probably 



of the most general interest since it contains the soft mode, is given in Table III. Total energy in mRy may be restored 
as the sum over all coefficients multiplied by corresponding powers of St in a.u. 

In constructing the kinetic-energy matrix and solving the eigenvalue problem for the TO T phonons in the harmonic 
approximation, we followed the guidelines of Ref. [Q (see also Ref. [Q for some details) . Calculated frequencies are 



given in Table [V along with the eigenvectors, which are expressed in terms of cartesian displacements of individual 
atoms (relative to the center of mass), multiplied by the square roots of masses, via back transformation from the 
symmetry coordinates defined in Table O. The comparison with experimentally measured frequencies is given in Table 
|y|. The overall agreement with experiment is within about 10%, with the exception of TO2 -Bi mode and the A2 mode. 
For the Bi phonon, the difference is probably due to some inaccuracy in the polynomial fit to the four-dimensional 
total energy hypersurface, since the frequency of the corresponding vibration (K - Nb stretching on the background 
of relatively undistorted oxygen octahedra) in the tetragonal phase (where it belongs to the Ai mode) came out in 
much better agreement with experiment in our previous calculations of Ref. P]. For the A2 phonon, the difference 
seems to be due to anharmonicity effects as has been already supposed in Ref. M. The explicit form of the total 
energy fit for this mode, which we found to be AE (mRy) = 15.16 (5*5)^ + 12.99 (S'5)^ (5*5 in a.u.), makes it possible 



to calculate the relevant frequency beyond the harmonic approximation, that is done in Section IV 



IV. QUANTUM OSCILLATOR TREATMENT 

Since the total energy as function of symmetrized displacements was available from our calculations in a sufficient 
number of points for providing the fourth-order polynomial fit, an extension of our analyzis beyond the harmonic 



approximation becomes feasible. Rather than introducing anharmonic corrections to frequencies obtained in the pre- 
vious section, we prefered to solve the Schrodinger equation for an, in principle, multidimensional quantum oscillator 
with a general-shape potential. We solved the Schrodinger equation numerically, by a finite-difference method (see, 
e.g., Ref. [|18|J19[ ), so that no basis functions have been employed. 

Similarly to how it was described elsewhere for the classical harmonic oscillator pq ] , we use an arbitrary convenient 
set of symmetry-adapted displacement coordinates St — X^i BuXi {xi are conventional cartesian displacements), which 
form a complete basis within a particular irreducible representation, but do not need to be orthonormal. Then, the 
Schrodinger equation aquires a form: 



|i:^a..4 + v((s.)) 



* = £;* 



where Gtt' — X^i Btim^ Btn is the kinetic-energy matrix. Mixed derivatives can further be excluded (that is desirable 
for maintaining high sparseness of the finite-difference matrix) by the following orthogonalizing transformation: 



-SSS*' 



Xt_ 
A 

where Xft is the t-th eigenvector, corresponding to the eigenvalue A*, of the kinetic energy matrix. 

The energy differences between the zero-point energy level and subsequent higher states give rise to characteristic 
frequencies of several lowest phonon modes. For the total energy fitted by the second order polynomial, these lowest 
frequencies should exactly coincide with the classical frequencies obtained in the harmonic approximation. This 
provides an additional mean of testing the accuracy of the finite-difference calculation scheme. Apart from allowing 
the anharmonicity in the most general way, the quantum treatment provides the description of high-frequency phonon 
branches (i.e. those lying beyond the lower 37V — 6 optical modes) and paves a way to incorporating the temperature 
effects. The energy levels of the oscillator may be further populated according to the Bose-Einstein statistics, and the 
transitions between levels analyzed. An analysis of this kind has been done by Bakker et al. [ pO[ for a one-dimension 
quantum oscillator model of the Ai phonon in LiTaOa. 

The accuracy of the calculation is essentially limited by the finiteness of the discrete mesh, and the convergence of 
eigenvalues in this parameter was explicitly controlled; however for the 4-dimensional oscillator {Ai and B2 irreducible 
representations), the maximal manageable grid of 14'*=38416 points was not sufficiently dense to produce as stable 
frequencies as those for Bi and A2 modes. 

For the B2 mode, the lowest energy levels of the 3-dimensional quantum oscillator were found to be 64, 89, and 130 
cm^^. Thus, the lowest possible phonon frequency is 25 cm^^, that is of correct order of magnitude, as compared 
with the experimental soft mode frequency (the harmonic approximation could only provide an information that the 
corresponding mode is indeed soft). Because of the small energy differences between several lowest oscillator levels, 
corresponding transitions can be easily excited already at medium temperatures, thus leading to a substantial increase 
of the soft mode frequency on heating. We however did not consider this effect in detail. 

For the one-dimensional A2 problem, the five lowest energy levels of the quantum oscillator are: 116, 357, 612, 880 
and 1160 cm^^. The lowest transition energy is therefore 241 cm~^ that is already an improvement over the result 
obtained in the harmonic approximation. The gradual increase of the separation between subsequent higher levels 
would result in somehow higher mean frequencies at elevated temperatures, but this effect is expected to be much 
smaller than for the soft mode. 

V. ATOMIC DISPLACEMENTS WITHIN THE SOFT MODE 

Now that the whole total energy hypersurface is mapped as function of symmetrized displacements of atoms 
from their equilibrium positions in the orthorhombic structure, it becomes possible to analyze the finite distortion 
compatible with the soft B2 mode which is a precursor of the orthorhombic to rhombohedral phase transition. We did 
not consider the change in the lattice strain which is an important clement of such transition. Nevertheless, it makes 
some sense to compare the symmetry-lowering displacements of atoms from equilibrium positions in the orthorhombic 
phase with the experimental geometry of the rhombohedral phase. 

We performed some additional fine-mesh fit near the position of the total energy minimum, in order to determine 
the ^-displacements of atoms with better accuracy than it is possible from the global fit by the 4th order polynomial. 
The equilibrium displacements of K, Nb and 0(1) in the [100] direction from the positions listed in Table | are 0.006 
A, —0.043 A and 0.013 A, correspondingly; they result in the total energy lowering of 0.06 mRy per formula unit. 



An additional degree of freedom enabling the displacement of Nb in the z direction along with that in x (i.e. the 
relaxation within both B2 and Ai modes) did not result in any better total energy minimum. The relative positions 
of thus displaced atoms projected on the (x,z) plane of the orthorhombic cell are shown in Fig. ^ in comparison with 
the experimentally determined positions in the rhombohedral phase. 

Magnitudes of the x-displacements of Nb and O in the rhombohedral phase are somehow larger than within the 
"frozen" soft mode. Otherwise, the most pronounced difference is the underestimated z-displacement of K, as has 
been already mentioned in section rt|. 

It is of considerable interest to know whether the displacements compatible with the B2 symmetry are energetically 
favourable if they occur only globally, simultaneously over the crystal, or also locally, on a short-range scale, as 
may be expected from the eight-site model (see, e.g., Ref. B for the discussion on the latter). We addressed this 



problem early for the case of tetragonal soft-mode displacements from the cubic phase |21|, and it was established 
that the [001] displacements result in a total energy lowering when they occur simultaneously in all unit cells, but 
not if they occur only in 25% of cells, as was tested in a supercell calculation. In other words, it is the indication 
that the cubic to tetragonal ferroelectric transition cannot be triggered by a spontaneous displacement of a single 
Nb atom, but rather needs a larger self-supporting region where such transformation occurs. [ P2[ Unfortunately, we 
could not give any estimates as for the size or shape of such regions from our supercell calculations. With respect to 
the orthorhombic phase of KNbOa, the similar theoretical study of the effect of local displacements was motivated 
by recent experimental observation by Dougherty et al. M| of relaxation processes in time-resolved Raman spectra 
which have been attributed to the hoppings, compatible with the Ai mode, over three classes of off-center Nb sites, 
distinguished by symmetry. There was however no experimental information available as for the spatial extent of 
spatial correlations in such relaxation processes. 

We tried to analyze several patterns of local displacements in the orthorhombic phase in supercell calculations, 
in order to establish some most important trends. Since it was not feasible to optimize the geometry in a supercell 
calculation in terms of all relevant atomic displacements, we limited ourselves to the study of the effect of pure 
Nb displacements which are obviously dominant. The incorporation of the relaxation of other atoms may somehow 
correct the magnitudes of the energy differences, but hardly the underlying qualitative trends. We considered local 
Nb displacements along z and x, i.e. those compatible with Ai and B2 modes, from the equilibrium position defined 
in Table ||. The relevant crystal superstructures modelled by 2xKNb03 and 4xKNb03 supercells are schematically 
shown in Figures and 0, along with corresponding total energy dependencies. 

Fig. |2| shows the effect on the total energy of z-displacements of every fourth "impurity" Nb atom, while other 
three Nb's in the supercell are kept fixed and provide a bulk macroscopic polarization. Apart from the minimum 
corresponding to Az(Nb)=0., the total energy grows rapidly; the displacement against the electrostatic field {Ai 
mode) through the centre of the Og nearest-neighbor octahedra to the symmetric position increases the total energy 
by about 30 niRy. Compared to the tiny value of 0.06 mRy energy lowering due to subsequent x-displacement {B2 
mode) into one of "eight sites" nearest to this point, it leaves no chance to find there a local minimum of the total 
energy with respect to both A^ and A^, variables. The contradiction with the experimental evidence of Ref. M 
that there are relaxation-type hoppings of the Ai symmetry between distinct local minima is apparently resolved 
by realizing that in the experiment, a considerable number of neighboring Nb atoms undergoes a relaxation hopping 
simultaneously, depending on the energy pumped in. A limiting case where really all atoms in crystal are involved 
in this process would mean the inversion of polarization, correspondingly mapping the total energy minima into new 
positions through the center of the unit cell. 

Another choice of 4xKNb03 supercell where "impurity" Nb atoms are organized in chains along the z direction (Fig. 
gb) leads to the same type of the energy dependence versus the displacement magnitude. It indicates that the local 
displacements against the bulk polarization are the main energetically unfavourable factors, whereas local frustrations 
within the z-chains of interchanging "impurity" and "bulk" atoms, which are peculiar for the first supercell geometry, 
are of minor importance. 

In contrary to this type of behavior, the x-displacements normal to the direction of macroscopic polarisation {B2 
mode) were found to be energetically favourable even if not supported by a similar displacement of near neighbors. 
This is in itself a clear indication of the order-disorder phase transition type from orthorhombic to rhombohedral 
structure. The precise energetics for different patterns of "up"- and "down"- Nb displacements from the (y,z) plane 
deserves however more detailed investigation. In addition to the perfect phase with the "frozen out" B2 phonon 
(Fig. ya), we considered two types of 2xKNb03 supercells which combine "up" and "down" displacements in two 
different ways: a chessboard arrangement with chains of identically shifted Nb atoms in the x direction (Fig. pp) 
and a NaCl-type structure where such chains are broken (Fig. yc). The corresponding total energy differences show 
that the last structure with frustrations in the x direction due to broken chains in absolutely unstable. Among two 
other structures which both contain undistorted chains, the one where "up" and " down" -chains appear in mixture is 
definitely more stable. An important consequence from this fact is its full agreement with the "chain structure" of the 



orthorhombic phase as proposed by Comes et al. |23|, whereas any other mixture of "up"- and "down" displacements 



which incorporates frustrations in the x direction can be relatively safely discarded. Another important concequence 
regards the mutual arrangement of chains. There seems to be a certain correlation within the (y,z) plane which favours 
each chain to have the neighboring chains with the opposite direction of displacement. In order to check that this 
statement is not confined to the particular chessboard geometry of the supercell in Fig. ^b, we performed an additional 
calculation where the effect of the displacement reversal within one chain out of four in the initial "perfect" geometry 
of Fig. ya has been studied. This was our possibly best supercell approximation to the displacement reversal within 
just one single chain in crystal at all, because thus reversed chains do not contact in this geometry. The result is that 
the spontaneous, say, "down" reversal of Nb displacements within one chain, with respect to the uniform background 
of "up" displacements, lowers the total energy by ~0.01 mRy per formula unit, i.e. is energetically favourable. The 
interaction between chains then works until some particular pattern with equal number of " up" and " down" chains is 
established. This stabilizes the orthorhombic structure, unless the lattice strain is allowed to become an independent 
parameter, and the system finds another global energy minimum compatible with the rhombohedral space group. 

Since the energy lowering due to x-displacements is smaller than the zero-point energy of the quantum oscillator 
in a potential well related to the B2 mode vibrations, it is worth noting that the above mentioned finite "up" 
and "down" displacements should not be understood literally, like well defined "frozen" atomic positions at zero 
temperature. Rather, we have to do with a damped quantum solution (in agreement with experimental observation 
of heavily damped B2 modes, Ref. [[7|) which is symmetric with respect to "up" and "down" x-directions for each 
single quantum oscillator. The interaction between chains distorts the symmetry of each individual solution due to 
the coupling between oscillators, and the mean value of corresponding x-displacement gets accordingly displaced. Its 
average over the whole system however remains zero. 

VI. SUMMARY 

We extended our previous theoretical study of TO F phonons in KNbOa to the orthorhombic phase. Based on 
the total energy calculations and polynomial fit of the results versus atomic displacements compatible with different 
irreducible representations, classical frozen phonon frequencies and corresponding eigenvectors have been calculated in 
the harmonic approximation, in satisfactory agreement with experimental data. For A2 and B2 modes, the quantum 
oscillator analysis has been done as well, resulting in some corrections to classical harmonic frequencies. The micro- 
scopic structure of the orthorhombic phase has been further analyzed in a series of supercell calculations modelling 
several patterns of Nb [100]-displacements compatible with the B2 mode. It was established that the formation of 
chains of uniformly displaced Nb atoms along the [100] direction is highly energetically favourable, in agreement with 
the experimental observations of the chain structure. The interaction between chains then favours the zero value of 
the average [100] Nb displacement, stabilizing the orthorhombic structure. 
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FIG. 1. Displacements of atoms, relative to the center of mass, in the (x,z) plane of orthorhombic KNbOa: from geometry 
optimization based on total energy calculations (black circles); from neutron diffraction measurements in the rhombohedral 
phase, Ref. pM (open circles). 

FIG. 2. Total energy differences (per formula unit) versus z-displacements of Nb from the equilibrium position for two 
superstructures. Thick line contours the supercell used. 

FIG. 3. Total energy differences (per formula unit) versus x-displacements of Nb from nominally equilibrium position in the 
orthorhombic structure for three superstructures. Thick line contours the supercell used. 



TABLE I. Positions of atoms in orthorhombic phase of KNbOa (in terms of lattice parameters) as determined by neutron 
diffraction measurements, Ref. [h4|, and optimized in the FP-LMTO total-energy calculation. 
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TABLE II. 

KNbOs. 



Symmetry coordinates used in the phonon calculations for four irreducible representations of orthorhombic 
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TABLE III. 4th order polynomial fit ol the total energy (in mRy) in terms ol the symmetry coordinates ol the B2 represen- 
tation (in a.u.). 
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TABLE IV. Calculated F-TO frequencies and eigenvectors in orthorhombic KNbOa 
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z 


0.08 


0.15 


0.25 


-0.37 


-0.37 


A2 


224 




X 








0.71 


-0.71 




146 


{ 


y 

z 


-0.85 


0.50 


0.00 


0.06 
0.08 


0.06 
-0.08 


-Bi 


232 


- 


y 

z 


0.07 


-0.09 


-0.63 


0.37 
0.39 


0.37 
-0.39 




297 


, 


y 

z 


-0.19 


-0.46 


0.67 


0.36 
0.16 


0.36 
-0.16 




528 


- 


y 

z 


0.13 


0.12 


0.26 


-0.37 
0.56 


-0.37 
-0.56 




82i 




X 


0.48 


-0.69 


0.43 


0.24 


0.24 


B2 


185 




X 


-0.73 


0.01 


0.32 


0.42 


0.42 




467 




X 


0.13 


-0.10 


-0.79 


0.42 


0.42 



TABLE V. Calculated and measured frequencies of F-TO phonons in orthorhombic KNbOa 



Mode 



Symm. 



Calculated 
frequency (cm~^) 



Expt; 



Expt.' 



Expt.'^ 



TOi 



TO2 



TO3 



TO4 



B2 
Bi 

B2 
Bi 

B2 
Bi 

Ai 

Bi 

A2 



^Raman spectroscopy; Ref. 

Raman spectroscopy; Ref. 

'^Infrared spectroscopy; Ref. 



soft 
232 
257 

185 
146 
186 

467 
528 
593 

297 
307 
224 



56 


40 


243 


249 


290 


281.5 


195 


196.5 


187 


192 


190 


193 


511 


513 


534 


534 


607 


606.5 


270 


270 


299 


297 


283 


283 



59 



197.5 



516 



A, (A) 



0.05 




K 



O 



0(1) 
0(11) 




K 



-0.1 -0.05 

# Nb 



Nb 




0.05 



0.1 A, (A) 



-0.05 








0.4 



0.8 A^ (a.u.) 




a 



b 



Fig. 2 



Ae 

(mRy) 




Ax (a.u.) 






a 



b 



Fig. 3 



